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ABSTRACT 

An  e3q>erimental  system  has  been  constructed  which  simulates  the 
detection  stage  of  a  receiver.  .With  the  application  of  periodic  bursts  of 
sine'wave  plus  -Gaussian  noise  at  the  input,  the  output  of  this  syster*  is 
non* stationary  random  process.  The  addition  of  noise  to  the  signal  siq>presses 
the  mean  at  the  output  for  the  half*wave  linear  detector.  However,  the  mean 
is  unchanged  for  the  half-wave  square-law  detector.  The  non -stationary 
variance  is  seen  to  vary  between  two  values  which  can  be  easily  measured. 
Output  signal-to-noise  ratio  is  greater  for  the  half-wave  square-law  device 
than  for  the  half-wave  linear  device.  Addition  of  a  clipper  to  the  half-wave 
linear  device  decreases  the  signal-to-noise  ratio.  Detection  characteristics 
of  the  half-wave  linear  and  half-wave  square-law  devices  are  found  to  be 
equivalent  and  the  clipper  degrades  the  detection  characteristic. 
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INTERACTION  OF  SIGNAL  AND  NOISE 
IN  A  NONLINEAR  SYSTEM 
by 

-P.  Trafton  and  D.  ' W.  Tufts 

Division  of  Engineering  and  Applied  Physics 
Harvard  University,  Cambridge,  Massachusetts 


I.  INTRODUCTION 

A  classical  problem  in  communication  theory  is  the  specification  of  the 
random  output  of  two  linear  filters  separated  by  a.  zero'-memory,  nonlinear 
element  when  excited  by  noise  and  signal  plus  noise,  as  shown  in  Fig.  1  below. 


s(t) +  n(t) 


FIG.  1.  Model  of  a  Nonlinear  Filter 


In  most  problems  of  interest,  the  first  filter  is  a  narrow-band  filter 
which  is  centered  at  a  high  frequency,  and  the  second  is  a  low-pass  filter. 

The  noise,  n(t)  ,  is  assumed  to  have  Gaussian  statistics,  and  the  signal,  s(t), 
is  often  a  burst  of  sine -wave. 

The  analytical  problems  encountered  in  this  type  of  problem  are  very 
formidable  and  only  limited  theoretical  results  are  available.  For  example,  the 
problem  of  finding  the  output  amplitude  distribution  function  is,  in  general,  un¬ 
solved.  In  most  cases,  the  first  two  moments  of  the  output  can  be  found  only 
approximately,  and  then  only  with  difficulty. 
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Considerable  theoiy  exists  so-called  envelope  detectors  with 
Gaussian  inputs.  In  this  situation,  the  output  of  the  narrow-band  filter  is 
represented  in  the  form*  E  ( t )  cos  ‘('a(t)  ]  »  where  is  the  center  fre- 

qmncy  of  the  filter.  The  random  envelope  and  random  phase  of  the  process  are 
E  ( t )  and  a  { t )  ,  respectively.  The  distribution  function  of  E  ( t )  is  known 
when  the  input  consists  of  Gaussian  noise  or  sine-wave  plus  Gaussian  noise. 

With  noise  alone,  E(t)  has  a  Rayleigh  distribution;  and  for  sine -wave  plus 
noise,  E  ( t )  has  a  modified  Rayleigh  distribution.  An  "envelope  detector "  is 
assumed  to  recover-  E  ( t )  or  some  power  of  E  ( t )  at  its  output.  Thus,  the 
output  statistics  are  known,  A  half-wave  linear  element  followed  by  an  ideal 
low-pass  filter  can  function  as  an  envelope  detector,  if  the  low-pass  filter 
passes  E(t)  exactly  and  rejects  all  harmonics  of  the  input.  This  theory,  how¬ 
ever,  is  not  applicable  if,  for  example,  the  low-pass  filter  does  some  filtering 
of  E  ( t )  .  Also,  this  theory  is  not  applicable  if  the  nonlinear  element  is  other 
than  linear  or  square-law.  Both  these  situations  are  explored  in  this  report. 

Power  spectrum  analyses  of  this  system  have  been  made  by  Rice  [1]  and 
Middleton  [8],  among  others,  and  the  essential  results  are  given  by  Davenport 
and  Root  [2].  These  analyses  lead  to  series  approximations,  in  most  cases, 
for  the  auto'-c  or  relation  functions  and  low-frequency  power  spectra.  From  this, 
the  first  two  moments  may  be  obtained.  Experimental  measurements  of  the  low 
frequency  auto- cor  relation  function  and  power  spectrum  for  a  variety  of  non- 
linearities  have  been  performed  by  Johnson  and  Fellows  under  the  direction  of 
Middleton  [8].  In  all  the  above  investigations,  the  input  is  often  Gaussian  noise 
or  sine -wave  plus  Gaussian  noise. 


TR450 


-3- 


Attempts  to  find  the  output  density  functions  have  been  made  in  a  series 
of  papers  by  Siegert,  and  Emerson  [3,  4,  5]  .  For  a  square >law  non> 

linearity,  Kac  and  Siegert  obtain  the  output  characteristic  function  but  inversion 
of  it  is  extremely  difficult.  Emerson  obtains  the  output  cumulants  and  'uses 
these  to  approximate  the  density  function  by  a  series  with  Gaussian  noise  and 
sine-wave  plus  Gat^ssian  noise  as  iiq>uts.  The  filters  are  assumed  t  have 
Gauss ian-shaped  amplitude  responses. 

In  view  of  the  analytical  difficulties  in  this  area,  an  experimental 
apparatus  has  been  constructed.  This  apparatus  permits  us  to  study  various 
attributes  of  the  random  output  of  Fig.  1  and,  in  particular,  the  first  two 
moments  and  the  amplitude  distribution  function.  A  wide  range  of  nonlinearities 
are  being  considered. 

n.  SYSTEM  DESCRIPTION 

The  diagram  in  Fig.  2  illustrates  the  e3q>erimental  system  that  is 
being  studied.  The  time  base  for  the  system  is  a  1  Kc  pulse  train  produced 
by  the  counter.  This  pulse  train  acts  as  an  external  trigger  for  the  pulse 
generator  which  then  produces  two  pulse  trains.  One  of  these  is  a  periodic 
1  Kc  sequence  of  rectangular  pulses  of  controllable  width  which  triggers  a 
gate  circuit.  Since  the  input  to  the  gate  is  a  continuous  130Kc  sine-wave, 
the  output,  s  (t) ,  of  the  gate  is  a  periodic  burst  of  sine-wave.  This  is  the 
signal.  The  second  pulse  train  produced  by  the  pulse  generator  is  delayed 
with  respect  to  the  first.  This  delayed  pulse  train  triggers  the  sampler  at 


V 


SYSTEM 
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the  output  of  the.Bystem.  Thus,  the  output  may  be  sampled  at  any  point  in 
time  with  respect  to  the  time  of  occurrence  of  the  leading  edge  of  a  signal 
pulse. 

The  pre-detection  filter  and  low-pass  filter  are  characterized  by  their 
amplitude  frequency  response  curves  and  examples  of  these  are  shown  in 
Figs.  3  and  4.  The  pre-detection  filter  consists  of  two  second-orde'r  tuned 
circuits  in  series  and  the  low -pass  filter  is  fourth-order  (SKL  Model  302). 
Both  filters  have  adjustable  bandwidths.  The  pre-detection  filter  performs 
an  adding  operation  so  that  its  output,  x(t) ,  is  a  sum  of  filtered  signal  plxis 
noise.  Also,  for  the  low -pass  filter,  |h(0)|=0;  i.  e.,  there  is  no 
transmission  at.  zero  frequency. 

The  stationary  noise  voltage,  n(t) ,  is  generated  by  a  gas  tube  and  is 
very  nearly  Gaussian.  It  is  not  exactly  Gaussian  because  extremely  large 
voltage  excursions  are  clipped,  an  unavoidable  effect  in  an  actual  physical 
system.  Its  power  spectrum  is  flat  across  the  frequency  range  of  the  pre¬ 
detection  filter. 

The  nonlinear  element  is  realized  by  a  so-called  photo-former  . 
This  is  a  unit  in  which  a  cathode-ray  oscilloscope  tube  (CRT)  is  the  main 
component.  The  following  diagram  illustrates  its  operation. 
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FIG.  3  FREQUENCY  RESPONSE  OF  LOW- PASS  FILTER. 
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FIG.  4  AMPLITUDE  RESPONSE  OF  PRE  -  DETECTION  FILTER. 
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FIG.  5  'Photo-former -Circuit 

The;CRT  spot  is  positioned  to.be  just  below  the  edge  of  the  mask.  If 
the.spot.moves  lower,  the  . increased  light  reaching  the  photo-tube  moves  the 
.spot  xtp.  If  the.  spot  moves  up  behind  the  mask,  the  reduced  light.moves  the 
spot. down.  The. spot,  therefore,  stays  at  the  edge  of  the  mask.  .As  the 
spot  is  .swept  horizontally  by  the  .oscilloscope  .sweep  circuit,  it  follows  the 
edge  of  the  mask,  tracing  out  the  pattern.  -Any  arbitrary  function  can  be 
generated  by  cutting,  out.  the  proper  mask  and  placing  it  over  the  CRT  face. 
The  input  signal  plus  .noise  is  then  applied  to  the  horizontal  input  and  the  out¬ 
put  waveform  is  then  obtained  from  the  vertical  output.  This  unit  has  been 
foiind  to  work  satisfactorily  if  the  input  frequency  is  moderate  (less  than 
200  Kc). 

The  elements  following  the  low-pass  filter  are  used  in  performing 
measurements  and  will  be  discussed  in  the  next,  section. 
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m.  MEASUREMENT  TECHNIQUES 

1.  Measurement  of  the  Amplitude  Diatribution  Function 

If  a.  fixed  periodic  eignal  plus  stationary  noise  is  applied  to  the  input 
of  the  system  of  Fig.  1,  then  the  output  is  a  non- stationary  random  process. 
Thus,  the  amplitude  distribution  function  varies  as  the  sampling  instant  1-^ 
varied,  and  in  our  case  (cf.  Fig.  2)  Uiis  variation  is  periodic.  The  measure' 
ment  is  performed  by  using  a  threshold  C(»nparator,  sampler,  and  counter  as 
shown  in  Fig.  2.  The  threshold  comparator  can  be  set  either  .to  pass  the  out> 
put  waveform  unchanged  except  for  an  arbitrary  dc  voltage  added  or  it 
si^plies  a  controllable  dc  voltage  to  the  sampler.  In  this  latter  case,  the 
sampler  (a  high-epeed  transistor  gate  of  2/isec  duration)  samples  a  dc 
voltage  at.  a.  1  Kc  rate.  The  threshold  on  the  counter  is  then  adjusted  to  cour.t 
at  this  level.  Then  the  random  waveform  is  .applied  to  the  sampler  and  the 
counter  indicates  the  number  of  samples  out  of  10,  000  that  exceed  the  thres' 
hold  (there  is  a  10  sec  counting  interval).  Thus,  we  have  an  experimental 
measurement  of  l''F(z)  at  some  time  t  from  which  F(z)  is  obtained  where 
F(  z)  is  the  amplitude  distribution  function.  By  varying  the  threshold  settings, 
we  obtain  F(z)  over  the  range  of  interest.  Then  the  sampling  tim^r  is 
changed  and  the  measurement  is  repeated.  Thus  we  measure  the  evolution 
in  time  of  the  amplitude  distribution  function. 

2.  Measurement  of  the  Mean 

In  the  presence  of  periodic  sigxial  plus  stationary  noise  at  the  input, 
the  mean  value  of  the  output  of  Fig.  1  also  varies  periodically  in  time.  The 
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measurement  of  the  mean  can  be  done  very,  simply  by  passing  the  samples  .of 
Fig.  2  directly  into  a  sensitive  galvanometer.  A  simple  calibration  is  performed 
by  passing  samples  of  fixed  and  known  amplitudes  into  the  galvanometer. 

3.  Measurement  of  the  Variance 

The  variance  of  the  output  of  Fig.  1  is  also  a  periodic  non> stationary 
function.  The  measurement  of  the  non>stationary  variamce  of  the  sa  aples  is 
somewhat  more  difficult  than  measurement  of  the  mean  or  amplitude  distribution 
function  and  requires  special  instrumentation.  This  ,  is  because  it  is  very 
difficult  to  find  an  instrument  that  can  perform  accurate  power  measurements 
on  pulses  that  are  very  narrow  compared  with  their  amplitudes.  Measure- 
ment  of  the  variance  is,  of  course,  a  power  measurement.  To  circumvent 
this  problem,  a  sample  and  hold  circuit,  (cf.  Fig.  2)  is  employed.  This  preserves 
the  amplitude  but  extends  the.  duration  of  the.  samples  and  it.  is  relatively  easy 
to  find  an  RMS  meter  that  can  perform  accurate  measurements  on  this  new 
waveform.  The  following  diagram  illustrates  the  operation  of  the  sample  and 


hold  circuit. 
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■Sampling  is  perfonned  once -every  1  m  sec.  The  sampling  duration 
is  2  and  tike  hold  duration  is  998 ji  sec.  Thus,  we  generate  pulses  of 

fixed  duratimi  and  random  amplitude. 

Ijet  us  -designate  the  value  of  z  (t )  at  each  sampling  instant  by 
where  i  designates  the  sampling  time.  This  value,  A.,  is  held  for  the 
-entire  period  of  the -output.  A.  is,  of  course,  a  randcnn  variable  ar-^  we  desire 
-to measure  its  variance)  i.  e. , 

N  N 

*•■*’'*  *  =  -r  Z 

i=l  i=l 

Following  well-known  methods  [1,  2],  we  can  easily  find  the  power 
--spectrum  of  the  output  of  the.  sample  .and  hold  circuit. 

p  (u)  =  ^  e  (u)  (3-2) 

where  T  s  1msec  in  our. system. 

S.iq>pose  now  that  we  put. this  waveform  into  an  RMS  meter.  What  does 
the  meter meaaure? -First  of  all,  the  delta  function  is  rejected  since  the  meter 
being  used  rejects  zero  frequency.  Thus,  the  meter  measures  the  following: 

The  integra.1  is  a  tabulated  one  and  equals  y  T/Z. 

A  RMS  =  =  \/var  (A)  (3-4) 
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Thus,  we  read  exactly  what  we  wish  since  var  (A)  is  the  variance  of  the 
sample  amplitudes.  Since  the  sample  and  hold  operation  may  be  done  At. any 
point  in  time,  the  variance  can  be  measured  as  a  function  of  time. 

The  accuracy  of  the  variance  circuit  can  be  easily  checked  by  the 
ergodic  theorem.  If  a  stationary,  zero  mean  noise  voltage  is  .applied  to  the 
RMS  meter,  the  time  average  of  the  variance  is  measured.  If  the  b'.me  noise 
voltage  is  passed  through  the  sample  and  hold  circuit.and  then  applied  to  the 
RMS  meter,  then  the  ensemble  average  of  the  variance  is  measured.  By 
ergodicity,  these  two  measurements  must  be  equal.  The  .sample  and  hold 
measurement  is  always  low  by  2  >3  %  .  Taking  this  into  account,  the  error 
in  the  variance  measurement  can  be  .made  very  small. 

It  may  be  noted  that  the  sample  and  hold  waveform  has  the  same 
statistics  as  the. samples.  Thus,  we  may  perform  further  measurements  on 
it  to  obtain,  for  example,  higher  moments. 

IV.  THEORETICAL  AND  EXPERIMENTAL  RESULTS  FOR  THE  MEAN 

It  is  possible  to  derive  some  analytical  results  concerning  the  behavior 
of  the  mean.  Consider  the  following  system:  * 


x(t) 


Zero -memory 
Nonlinear 
^  Element 


>y(t) 


FIG.  7  Nonlinear  Element 
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Astvming  .x(t)  to  be^the  ou^ut  of  a  narrow-band  filter  centered  at 
fraqnenc^r  f^  ,  we  can  eaq»reee  xCt)  as  follows: 

x(t)  s  E(t)cos  [2vf^t.4a  (t)]  (4-1) 

E(.t)  is  tbe.random  envelc^  and  a  (t)  the  random  phase  of  x(t)  .  We 
.r^resent  the  nonlinear  operation  ns  .follows: 

y  =  f(x)  (4-2) 


or 


j  i.(Ecos0.')  where  B  -  2w  f  t-+ a  (4-3) 

o 

.FoUowini;  a.  currently,  popular  techniq^  [6,7]  ,  we  eaqtand  y  into  a  . Fourier 
cosihe.series  since  y  .is  an  even  periodic  function. of  0. 


y 


1 

T 


a_  +,a.  COS  0  + 


a^cos 


20  + 


and 


(4-4) 


f  (£  cos  0  )  cos  n0  d0 


(4-5) 


/We. are  interested  in  the  a^  term  in  y  .since  the  ou^ut. is, terminated  by  a 
low-pass  filter. 

Consider  now  a- specific  nonlinearity,  namely,  a  half-wave  linear 
.element.  This  device. is  specified  as  follows; 

f  (E  cos  0  )  =  E  cos  0  ,  E  cos.0  >  0 

(4-6) 

f  (£  COS.0  )  -  0,  E  cos  0  <  0 
and,  of  course,  E(t)  >  0  . 
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*o  = 

4.r  • 
0 

cos  0  60  ~ 

(4-7) 

«y(t) 

- 

(Ignoring  high-frequency  terms) 

(4-8) 

y(t) 

.  E(t) 

» 

(4-9) 

where  .the  overb&r  indicates  statistical  averaging. 

If  the  input, noise  to  the  i»rr6w«band  filt&r  is  Gaussian,  then  £(t) 

,is  known  to  have  a  Rayleigh  distribution  in  the  presence  of  noise  only.  In  the 
presence  of  sine«wave  plus  -Gaussian  noise,  E  ( t-)  is  known  to  have  a  modified 
Rayleigh  distribution.  The  means  for  both  aituations  are  known. 

In  the  absence  of.sine-wave: 


E{t)  = 


(4.10) 


In  the  presence  of  sine^wave  : 


E(t)  = 


(re 


where  A  is  the  amplitude  .of  the  input  sine  .wave  and  v  is  the  RMS  noise 
power  at.  the  input  of  the  nonlinearity.  I^  and  I^  are  modified  Bessel 
functions  of  orders  zero  and  one,  respectively 
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At  this  point,  it  is  convenient,  to  define  .a  signal-to-uoise- ratio.  Let 
p  -  .  This  gives  the  ratio  of  the  .mean  square  of  the  sine-wave  to  the 

jnesA  square  of  the  noise.  Therefore, 

_  -4 

W*>  =  *  t(l+P  )  lo*  -f  >  +  P  <  ■§•  )  1  <4-14) 

y£g 

In  our  experimental. apparatus,  the  time  interval,  during  which  the  sine- 

; 

wave  is  on  is  small  compared  to  the  period  of  the. signal  (100 /i  sec  vs.  1msec). 
On  a  time  average. basis,  therefore,  the  mean  is  very  nearly  y^(t)  . 

However,  the  low-pass  filter  does  not  transmit  zero  frequency, and_so  the  mean 
at  the  filter  output,  is  . the  response  of  the.filter  to  the  .difference  of  the  mean 
with  and  without  signal..  The  following  diagrams  illustrate  this  behavior. 


FIG.  8  Meanofy(t)  Ignoring  High-Frequency  Terms 
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Subtracting  out.  flie  .  lime  average  of  the  mean,  we  ojbtain: 


L/' 


FIG.  9  UTaveform  of  FIG.  8  Minus  y  (t) 

Let.  the  output  of  the  low-pass  filter  be  .z  ( t )  .  Therefore, 

00 

*(t)  =  J  y(T)h(t-T)dT  (4-15) 

-oo 


and 


00 

*(t)  =  J  y(T)  h(t-T')  dr 
-oo 


where 


(4-16) 


y(T)  =  Yg+n^r)  -  Y^ir)  (4-17) 

An  interesting  behavior  of  the  mean  for  a  half-wave  linear  ,  element 
can  now  be  demonstrated. 

_  -4 

y.(T)  ^  -sr*  ^  [(1+P)  I  ( 4>  +P  I.  (4-)]  -  — 

\/2t  /Z? 


(4-18) 
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-Siq>po8e  now  that  p~*’ao,  i.  e. ,  o*  0  .  Using  the  asymptotic :behavior 
of  the  modified  Besael  functions  and  ,  it.  can  be  easily  shown  that 
y  (t).’*’  fids  pulse  is  the  input  mean  to  the  low-pass  filter,  then  we 

call  the  output  mean  the  mean  without,  noise.  Further,  yCr  )  is  .monptonic 
in  p  .  Therefore,  if  we  fix  A  and.s  is  either  zero  or  some  non-zero 
constant,  then  yTrT  with  noise  (s  /  0)  is  less  than  yCr)  .without -oist 
(o’  -  0).  This  is  a  sufficient  condition  to  prove  that  the  mean,  z  (t)  ,  at  the 
output.of.the  low-pass. filter  is  always. suppressed  in  the  presence  of  noise. 

U  .z^(t)  .is  the  output  mean  in  the  absence  of  noise  and 
put  mean  in-the;presence  of  noise,  then  the  following  is.  true: 

l^^ll  <  l*,(t)|  (4-19) 

and  furthermore,  ~  where  K  is  less  than  one  and  can  be 

-calculated.as  a  function  of  p  . 

InJFlgs.  14a,.b.at  the.  end  of  this  report,  experimental  curves  are  shown 
demonstrating  this  effect  along  with  experimental  and  calculated  values  for  K, 
the  aiq>pres8ion  factor.  The  low-pass  filter  -has  been  set  to  give  considerable 
filtering  of  the. signal  envelope.  P.re-detection  3  dB  bandwidth  is  15 Kc  and 
low-pass.  3  dB  bandwidth  is  8  Kc. 

Let  us  now  cons^er  the  .case  .where  the  nonlinear  element  is  a  half- 
wave,  square -law  device.  In  this. case, 

f  (E  cosd)  =  cos^^  ,  E  cos  0  >  0 

(4-20) 

f  (E  cos  6  )  =  0  ,  E  cos  0  <  0 
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Proceeding  as  before,  we  find 


*io  '  2  * 


(4-^21) 


1  ^2 


y{t)  =  -^.E  (t)  (Ignoring  high-frequency  terms) 


In  the  absence  of  sine -wave; 


- r:^  1  -  2.  ,  a 

y  ( t)  =  -r  ( t)  =  — :t- 
'n'  4  n  '  '  2 


In  the  presence  of  sine-wave  : 


(4-22) 


(4-23) 


(4-24) 


(4-25) 


This  leads  to  an  interesting  resvdt.  since  the  above  is  independent  of  the  noise. 
Thus,  the  output  mean  is  the  same  both  with  and  without  noise,  i.  a. , 

*.+n<‘>  =  *8<‘> 

In  Fig.  15  at  the  end  of  this  report,  experimental  curves  are  hown 
demonstrating  the  above  result. 

Now  suppose  that  a  half-wave  linear  with  clipping  characte  r Isuc  is 
used  as  follows:  y 


FIG.  10  Half-wave  Linear  Characteristic  with  Clipping 
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Proceeding  as  before,  we  find 

r(t)  =  -f-  E(t),  0<  E(t)<  N 

y(‘)  =  ,  E(t)  >  N 

The  density  functions  of  £  ( t )  with  and  without  signal  appear  as 

follows 


FIG.  11  Density 'Functions  of  £  (t).  With  and  Without  Signal 

Inspection  of  this  diagram  shows  that  ~  is  less  now  than  when 

the  clipper  is  absent  although  the  amount  cannot  be  calculated  exactly.  This 
indicates  additional  suppression  of  the  mean  with  the  addition  of  the  clipper. 

InFigs.L4a,b  data  showing  this  effect  are  given.  The  clipper  passes 
signal  without  noise  xinchanged.  Therefore,  the  curves  of  mean  without  noise 
are  identical  to  those  for  the  half-wave  linear  device.  The  mean  with  noise, 
however,  is  somewhat  further  suppressed  below  that  of  the  vmclipped 
detector. 

It  should  be  pointed  out  that  the  dc  rejection  of  the  low-pass  filter 
does  not  produce  the  above  effects  but  instead  displays  them.  The  important 
point  is  that  the  time-varying  response  of  the  mean  behaves  as  shown  regardles 
of  any  dc  values  that  may  or  may  not  be  present. 
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V.  THEORETICAL  AND  EXPERIMENTAL  RESUL'K  FOR  THE  VARIANCE 

In  the  presence  of  a  burst  of  sine-wave,  the  ou^ut  of  the  nonlinear 

element,  y(t),  (cf.  Fig.  2)  has  a  non-stationary  auto-correlation  function 

which  we  .designate  as  R  (t.  ,  t.)  .  However,  if  we  concern  ourselves  with 

y  - 1  c 

the  low-frequency  part  of  y(t),  namely  '*£(t  )  or  Bome  power  or  modified 
version  of  it,  then  E  ( t)  is  a  stationary  process  for  noise  alone  as  input  or 
for  sine-wave  plus  noise. as  input.  Thus,  the  low-frequency  output  is,  so  to 
speak,  "piecewise”  stationary^  The  following  diagram  of  the  t^  ,  t^  plane 
.shows  the  behavior  of  the  low-frequency  a.uto-correlation  function  as  we  switch 
from  noise  to  signal  plus  noise. 


FIG.  12  Diagram  Pertaining  to  a  Comparison  of  R  (t,  ,  t_  ),  With  and.  Without 
-Signal  y  ^  ^ 
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Sine-wave  is  turned  on  at  t^  =  t^  =  0  .  designates  that 

s(t2)  =  0,or  vice-versa.  This  behavior  occurs  in  the  first  and 
third  quadrants.  R^Ct^jt^)  is  stationary  in  each  quadrant  of  the  t^,  t^ 
plane,  i.  e. ,  R^ltj  ,  t^)  =  Rytr)  where  ~  *2  *  *1  "  express  the 

variance  at  the  output  of  the  low-pass  filter  as  follows: 

oo  oo 

var  -tj)h(t-t2)dtj  dt^  (5-1) 

-oo  -oo 

In  problems  of  this  sort,  it  can  be  shown  that  the  auto-correlation 
function  can  be  written  in  a  series  [2],  the  first  term  of  which  produces  the 
squared  mean  at  the  low-pass  ou^ut.  In  what  follows,  we  assume  that  this 
term  has  been  subtracted  out,  so  that  the  above  gives  the  variance  and  not 
the  mean  square  of  the  output. 

We  can  rewrite  £q.  (5-1)  as  follows: 

OD  oo  oo  oo 

var  ( t  >]  =  I  -t^)h{t-i2)dtjdt2 

-00  -00  o  o 


oo  o 


-2 


I  I  R'g+„(T)h(t-tj)h(t-ydtjdt2 


(5-2) 


O  -OO 

The  iirst  term  may  be  rewritten  as 
oo 

y  P_j«)  |h(w)|^  du 
-oo 


(5-3) 
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wliere  P  (u)  and  H(y)  are  the  Fcmrier  tranafonna  of  R  (r)  and  h.{r)  , 
n  n 

re^ectieely.  Thia  term  ia  recognised  aa  tiie.  ateady-atate  reaponae  of  die 
variance  at  the  loar>paaa  ott^ot  when  only  noiae  ia  preaent  in  the.  ayatem. 
The  aecond  term  may  be  written  aa 


QD 


(5-4) 


where 


=  j"  h{x)e*j^  ix 


(5-5) 


Thia  term  ia  a  fonction  of  time  aa  ia  the  third  term. 

For  t  <  0  ,  the  variance  ia  given  by  the  following 


oo 

{»(t)]  =  j  P^(u)|H(u)|^<iu 


(5-6) 


Thia  ia  the  atationary  variance  of  the  n( t)  proceaa  at  the  output  of  the  low> 
paaa  filter .aaauming  that  n(t)  haa  exiated  at.the,iiq>ut  for  a  time  interval 
greater,  than  the  memory,  of  the  low-paaa  filter. 

Aa  t  oo  ,  aaauming  that  the  a(t}  +  n^t)  proceaa  ia  left  on  for 
a  time  greater  than  the  memory  of  the  low-paaa  filter,  the  variance  is  given 
by  the  following 


oo 

=  y  1h(w)|^  du 


(5-7) 
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This  is  the  stationery  variance  of  the  s(t)  +  n(t)  process  at  the  output  of 
the  low-pass  filter. 

The  low-frequency  variance  at  the  output  of  the  half-wave  nonlinear 
olement  will,  in  general,  have  the  following  form 


FIG.  13  Variance  of  y(t)  Ignoring  High-Frequency  Terms 

If  the  low-pass  filter  had  sufficiently  wide  bandwidth,  then  the  output 
variance  would  have  the  same  form  as  above,  i.  e. ,  the  low -pass  filter  would 
completely  pass  the  low-frequency  spectrum.  However,  if  the  low-frequency 
spectrum  is  filtered,  then  this  form  will  not  he  preserved.  Experimental 
data  for  the  variance  are.ahownin  Fig.  16  for  the  half-wave  linear  element. 

If  only  noise  is  present  at  the  input,  then  the  output  variance  is  given  by 
£q.  (5-6)  and  is  shown  experimentally  by  the  line  labeled  "noise  alone.  " 

If  the  signal  were  left  on  indefinitely,  i.  e. ,  the  input  is.  sine-wave  plus  noise, 
then  the  output  variance  wovild  be  a  constant  given  by'Eq.  (5-7)  and  is  shown 
experimentally  by  the  line  labeled  "sine -wave  plus  noise.  "  Thus,  the  non¬ 
stationary  variance  is  seen  to  equal  the  value  given  by  Eq.  (5-6)  before  the 
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burst,  occurs  .and  then  increases  toward  the  value  given  by;  £q.  (5-7) 
when  the.signal  Jburst  is  applied.  Then,  when  the  signal  is  turned  off,  the 
variance  returns  to  the  value  given  by:£q.  (5-6). 

If  the.  signal  burst,  were  left  on  for  a  aufficient.  length  of  time,  then  we 
would  e^ect  the  non- stationary  variance  to  attain  the  value  given  by.  £q.  (5-7). 
That  this  actually  occurs  is  shown  in  Fig.  17  where  the  burst,  was  left  on  for 
a  longer  time  but.  all  other  parameters  were  .unchanged. 

In  Figs.  18  and  19  are^shown  the  variance  curves  for. the  half-wave 
^.linear  with  clipping  and  half-wave -square  law  devices.  The  .clipper  reduces 
the  variance,  as  expected,  and  the  half-wave  square  law  device  exhibits  a 

very  pronounced  increase  in  the  variance  during,  the  time- varying  response. 

* 

VI.  TH£  M£AN  SQUARE 

A  measurement  of  the  mean  square  has  not  been  considered  necessary. 
The  .mean  sq\iare  of  the  output,  z  ( t ),  is  given  by  the  following: 

A  X 

z  (t)  =  .z(t)^  +  var{z(t)]  (6-1) 

.  Since  the  mean  and  the  variance  .can  be  accurately  measured,  the  mean  square 
may  be  calculated.  Experimental  data  on  the  mean  square  are  given  in 
Figs.  20,-  21, .  22.  Since  the  squared  mean  is  positive  and.  the  variance. was 
not  observed  to  go  below  its  value  for  noise  alone  as  input,  it  is  seen  that 
the  output  mean  square  is  increased  in  the  presence  of  signal.  This  can  be 
used. as  a  basis,  for  detecting  a-signal. 
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Vn.  AMPJUTUDE  DISTRIBUTION  FUNCTION 

Purely  . theoretical  results  for  the  amplitude  distribution  function  are 
difficult. to  obtain.  This  is  the  problem  of  passing  a  non-Gaussian  process 
throu^  a  linear  filter  and.calculating  the  output  distribution  function.  This 
is  an  unsolved  problem. 

At.  this  -stage  of  our  .-work,  the  amplitude  distribution  functions  have 
been  measured,  for  the  nonlinearities  already  mentioned.  The  .envelope  has  a 
Rayleigh  density  but. with  filtering  of  the.  low-frequencies,  the  density  is 
unknown. 

Figures  23a,  23b  and  24.show  the  measured  amplitude  distribution 
functions,  or,  to  be  exact,  <me  minus  .the  distribution  function.  At  180  p  sec 
delay. with  respect  to  the  leading  edge  of  the  input  pulse,  the. effect  of  the 
.signal,  is  .not  present.and ..the  distribution  function  of  noise  is  measured.  At 
280 It  sec  delay,  the  maximum  of  the  output. mean  occurs.  The  distribution 
function  is.shifted.to  the  right  due  to  the. change  in  the  mean  and  also  is  seen 
to.have  an  .increased. spread,  indicating  an  increase  in  the  variance.  Both 
of  these  .observations  are  consistent  with  the  data  for  the  mean  and  variance. 

If  Gaussian  noise  is  passed  through  the  system  and  a  half-wave  linear 
element,  is  used,  then  we  should  obtain  a  Rayleigh  density  at  the  output  if 
the  low-pass  bandwidth  is  sufficiently  wide.  If  a  half-wave  aquare-la.w  device 
is  used,  then  an  exponential  density  is  expected.  These  distributions  have 
been  obtained. and  are  shown  in  Figs.  25  and  26.  Also  shown  here  are  plots 
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of  tot-fa|_l-F(*)j,]  /VS.  biz  .  For  a  ^yleigh  distributioiiy  tliis  vrill  give  .a 
■  Straight  line  .wiih  slope  2. 00.  For  an  eiq»onential  distribution,  it.  will  give 
.a:straight.line.witb.slope,l.  00.  These  data  dem<mstrate  that  the.system  can 
produce  bnown  results. 

Vm.  OUTPUT  SIGNAL-NOISE  RATIO  AND  DETECTION  CHARACTERISTICS 
‘Signal-noise  ratio  is. not. a  unique  .quantity  but  the  fcdlowing.  choice  .has 
been.studied.in  this  work  *. 


SNR  s 


n 


(8-1) 


is  the  .output; when  signal  plus  noise  is. applied  at  the  input  and 
Z]^(  t)  is.the  output  when  noise.alone  is  present.in  the.system.  This  .choice 
of  ■  SNR  thus  gives  a  measure  of  the. excess  power  in  the  output.when  signal 
.is  applied,  to  the  .input.  Curves  in  Figs.  27,  28,  and  29.  show  the  output  SNR 
for  three  nonlinearities.  The  balf-wave .square-law  device  gives  a  .higher 
■  SNR  than  the  .half-wave  linear,  and  the  presence  of  a  clipper  on  the  half-wave 
linear  degrades  the  .SNR. 

From  a  detection  point  of  view,.  SNR  is  not  necessarily,  significant 
and  the  detection  characteristic  is  ol‘  greater  interest.  This  gives  the 
probability  of  detection  for  a  given  false  nlarm  p  robability  and  this  data  can 
be  read  directly,  from  the  amplitude  distribution  curves.  Detection 
characteristics  .are  given  in.Figs.  30  and  31.  These  .ci  rves  show  that 


TR450 


-24- 


the  half-wave  square  law  and  the  half-wave  linear  devices  are  equivalent,  and 
the  half-wave  linear  device  with  clipping  degrades  the  detection  characteristic. 
This  cannot  be  shown  analytically,  since  the  output  distribution  functions  cannot 
be  calculated.  The  above  can  be  demonstrated  analytically  for  strictly  envelope 
detectors,  however,  for  detecting  signals  of  tinknown  phase  [9]. 

DC-  FUTURE  WORK 

The  experimental  system  at  present  simulates  detection  plus  low-pass 
integration.  Other  forms  of  integration  will  be  examined,  in  particular,  RC 
integration  and  ideal  integration,  i.  e. ,  unweighted  integration.  The  effects 
of  these  forms  of  integration  on  signal-to-noise  ratio  and  signal  detectibility 
will  be  examined.  The  effect  of  different  nonlinear  elements  will  also  be 
investigated. 

An  attempt  to  classify  the  distribution  functions  at  the  output  of  the 
low-pass  Rlter  and  integrator  will  be  made  using  the  Pearson  class  of 
distributions  and  also  series  approximations  in  terms  of  moments.  To  this 
end,  the  third  and  fourth  moments  will  be  measured. 
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MEAN  VS.  TIME 


FIG.  14a  MEAN  RESPONSE  FOR  HALF-WAVE  LINEAR  DEVICE 
WITH  AND  WITHOUT  CLIPPING  WITH  CONSTANT 
NOISE  POWER  (a^) 


SAMPLE  MEAN  (VOLTS)  SAMPLE  MEAN  (VOLTS) 


FIG.  14b  MEAN  RESPONSE  FOR  HALF-WAVE  LINEAR  DEVICE 
WITH  AND  WITHOUT  CLIPPING  WITH  CONSTANT 
NOISE  POWER  (0*2) 
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MEAN  vs.  TIME  (WITH  AND  WITHOUT  NOISE) 
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FIG.  15  MEAN  FOR  HALF-WAVE  SQUARE  LAW  DEVICE  WITH 
INPUT  NOISE  POWER  CONSTANT 
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WITH  INPUT  NOISE  POWER  CONSTANT  (ff2) 


RMS  OEVIATiON  FROM  MEAN  VS.  TIME  (  PULSE  DURATION  -  ISO  |1SEC ) 

SMRv98  _ 


FIG.  17  VARIANCE  WITH  EXTENDED  PULSE  DURATION  AND 
CONSTANT  NOISE  POWER  (ff^) 
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WITH  INPUT  NOISE  POWER  CONSTANT  (0^) 
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FIG.  19  SQUARE  ROOT  OF  VARIANCE  FOR  HALF-WAVE  SQUARE  LAW 
DEVICE  WITH  INPUT  NOISE  POWER  CONSTANT  (O^) 
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CONSTANT  INPUT  NOISE  POWER  (02) 


WITH  CONSTANT  INPUT  NOISE  POWER 
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FIG.  22  MEAN  SQUARE  FOR  HALF-WAVE  SQUARE  LAW  DEVICE 
WITH  CONSTANT  INPUT  NOISE  POWER  (O^) 


FIG.  23a  AMPLITUDE  DISTRIBUTION  FUNCTION  FOR 
HALF-WAVE  LINEAR  DEVICE. 
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FIG.  24  AMPLITUDE  DISTRIBUTION  FUNCTION  FOR 
HALF-WAVE  SQUARE  LAW  DEVICE. 


PRE- DETECTION  BANDWIDTH  •  S  Kc 
LAWS  PASS  BANDWIDTH  •  12  KC 


OO 


Q 


FIG.  25  EXPERIMENTAL  RALEIGH  DISTRIBUTION 


FIG.  26  EXPERIMENTAL  EXPONENTIAL  DISTRIBUTION 
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FIG.  28  OUTPUT  SNR  FOR  HALF-WAVE  LINEAR  DEVICE 
WITH  CLIPPING. 
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FIG.  30  DETECTION  CHARACTERISTICS  FOR  HALF-WAVE  LINEAR 
ANO  SQUARE  LAW  DEVICES. 


FIG.  31  DETECTION  CHARACTERISTICS  FOR  HALF-WAVE  LINEAR 
DEVICE  WITH  CLIPPING. 
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